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Event-by-event fluctuations of the chemical composition of the hadronic final state of relativistic 
heavy-ion collisions carry valuable information on the properties of strongly interacting matter 
produced in the collisions. However, in experiments incomplete particle identification distorts the 
observed fluctuation signals. The effect is quantitatively studied and a new technique for measuring 
chemical fluctuations, the identity method, is proposed. The method fully eliminates the effect of 
incomplete particle identification. The application of the identity method to experimental data is 
explained. 

PACS numbers: 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

Event-by-event fluctuations of chemical (particle-type) 
composition of hadronic final states of relativistic heavy 
ion collisions are expected to be sensitive to properties of 
strongly interacting matter produced in the collisions [l[ . 
Specific fluctuations can signal the onset of deconfinc- 
ment when the collision energy becomes sufficiently high 
to create droplets of quark-gluon plasma |2(. At higher 
collision energies, where the quark-gluon phase is abun- 
dantly produced at the early collision stage, large chem- 
ical fluctuations can occur as the system hits the critical 
point of strongly interacting matter in the course of its 
temporal evolution [1, 0|- It is thus certainly of inter- 
est to study event-by-event chemical fluctuations exper- 
imentally. First data coming from the CERN SPS 0-0] 
and BNL RHIC [1] were already published. The results 
are not very conclusive yet and more systematic mea- 
surements are needed. In addition the question arises 
whether data analysis methods can be improved. 

In real experiments it is impossible to determine 
uniquely the type of every detected particle. The identi- 
fication requires measurements of particle electric charge 
and mass. Precise mass measurements are experimen- 
tally difficult and expensive. For this reason analyzes of 
chemical fluctuations are usually performed in a limited 
acceptance where particle identification is relatively reli- 
able. However, sensitivity to fluctuations of range larger 
than the acceptance window is then lost and signals from 
fluctuations of shorter range are usually diluted. Further- 
more it should be noted that results on fluctuations, un- 



like those on single particle spectra, cannot be corrected 
for the limited acceptance. Often it is possible to enlarge 
the acceptance, but at the expense of a significant con- 
tamination of the sample by misidentified particles. The 
effect of particle misidentification can distort measured 
fluctuations. Thus, incomplete particle identification is a 
serious obstacle to the precise measurement of chemical 
fluctuations. 

Although it is usually impossible to identify every de- 
tected particle, one can in general determine with high 
accuracy the percentage (averaged over many interac- 
tions) of, say, kaons among produced hadrons. This in- 
formation will be shown to be sufficient to fully eliminate 
the effect of incomplete identification. In this paper we 
propose a new experimental technique called the identity 
method which achieves the goal independently of the spe- 
cific properties of the chemical fluctuations under study. 

In the study of chemical fluctuations the NA49 Col- 
laboration [5H7( used the measure Odyn, which is defined 
as the difference between fluctuations measured in real 
and mixed events. The effect of particle misidentifica- 
tion is accounted for by including it in the mixed events. 
The STAR Collaboration Q used, in addition to the (Tdyn 
measure, the quantity t'dyn- The latter one assumes that 
particles are uniquely identified. 

It was suggested long ago [j| [Tfl to quantify chemical 
fluctuations by the measure $ [ll[ which proved to be 
efficient in experimental studies of event- by- event fluctu- 
ations of particle transverse momentum [l3 . [l3j , electric 
charge and quite recently of azimuthal angle fl5j . 
The <f> measure, unlike cxdyn and v^yn, is a strongly inten- 
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sive measure of fluctuations. Namely, its magnitude is 
independent of the number and of the distribution (fluc- 
tuation) of the number of particle sources, if the sources 
are identical and independent from each other. This fea- 
ture, which is discussed in detail in [16j . is important in 
experimental studies of relativistic heavy-ion collisions 
where the collision centrality is never fully controllable. 
However, up to now it was unclear how to correct mea- 
surements of chemical $ for the effect of particle misiden- 
tification. 

The identity method, which is developed here, uses the 
fluctuation measure "J, a simple modification of $. The 
measure "J, similarly to $, is strongly intensive but the 
modification allows to correct the measurements for the 
effect of particle misidentification. Below we show that 
the measure can be factorized into a coefficient, which 
represents the effect of misidentification, and the quan- 
tity "J ci which corresponds to the value VP would have for 
complete identification. The misidentification coefficient 
can be determined from the data in a model independent 
way. Therefore, the identity method provides the value 
of the fluctuation measure as it would be obtained in an 
experiment in which every particle is uniquely identified. 

Before the identity method is presented we introduce 
and discuss in Sec. |H] the fluctuation measures Odyn, ^dyn, 
<I>, and In Sec. IIIII we demonstrate by a Monte Carlo 
simulation how the effect of misidentification distorts the 
chemical fluctuations as quantified by i/^yn, $ and 
The identity method is formulated in Sec. IIVI Instead 
of conclusions we present in the last section the steps 
required to apply the identity method to experimental 
data. In order to simplify the presentation, we consider 
chemical fluctuations of events composed only of kaons 
and pions. Clearly, kaons and pions can be replaced by 
particles of any other sort. 

II. MEASURES OF FLUCTUATIONS 

As mentioned in the Introduction, fluctuations of 
chemical composition of final states of relativistic heavy- 
ion collisions can be studied in several ways. The NA49 
Collaboration [5 7] measured event-by-event fluctuations 
of the particle ratios K/ir, K/p, p/ir and determined the 
quantity <7d yn defined as 



CTdyn = Sgn((Td ata - CTmixcd) Y kdata ~ ^mixcdl > W 

where <7d a ta and <7 m i xe d is the relative width (the width 
divided by the mean) of the event- by- event particle ra- 
tio distribution in, respectively, the data and artificially 
generated mixed events where every particle comes from 
a different real event. The fluctuations present in mixed 
events are due to the effect of particle misidentification 
and the statistical noise caused by the finite number of 
particles. 

The STAR Collaboration used [|| the quantity v^yn to 
measure chemical fluctuations. For the case of a two- 



component system of pions and kaons Vdyn is defined as 

(N K (N K - 1)) {NAN, - 1)) _ (N K N V ) 
Udyn (N K ) 2 (N^y (N K )(N„) ' 

(2) 

where Nx and N v are the numbers of kaons and pions in 
a given event and (...) denotes averaging over events. The 
quantity z^dyn is defined in such a way that, in particular, 
it vanishes when the multiplicity distributions of pions 
and kaons are both poissonian ((Ni(Ni — 1)) = (A^) 2 , 
i = 7T, K) and independent from each other ((NkN,) — 
(Nk)(N„)). Thus, it is constructed to quantify the de- 
viations of the fluctuations from the poissonian noise. 
For large enough particle multiplicities one finds the ap- 
proximate relation v dyn » errata ~~ '''mixed which gives 
j'dyn ~ sgn(iTdyn) o"dyn H ■ We note that the quantity 
fdyn implicitly assumes unique identification of all parti- 
cles. 

As already noted, it was advocated long ago [{| 
to employ the measure $ [ll| to study chemical fluctua- 
tions. The measure is defined in the following way. One 
introduces the variable z = x — x, where x is a single 
particle characteristic such as the transverse momentum 
or azimuthal angle. The over-line denotes averaging over 
the single particle inclusive distribution. The event vari- 
able Z, which is a multi-particle analog of z, is defined 
as Z = Y^i( x i ~ ^Oi where the sum runs over the N 
particles in a given event. By construction, (Z) = 0. The 
measure $ is finally defined as 

The measure $ vanishes in the absence of inter-particle 
correlations. This situation is discussed in some detail 
below for the case of chemical fluctuations. Here we note 
that the measure also possesses another important prop- 
erty - it is strongly intensive which means that it is inde- 
pendent of the number and of the distribution (fluctua- 
tion) of the number of particle sources, if the sources are 
identical and independent from each other. In particu- 
lar, if a nucleus-nucleus collision is a simple superposition 
of nucleon-nucleon interactions, then $,4,4 = <5>nn- The 
strongly intensive property is a very valuable feature of 
$, as centrality selection in relativistic heavy-ion colli- 
sions is never perfect and events of different numbers of 
particle sources are always mixed up. The strongly inten- 
sive property is also desirable when different ccntralities 
or different colliding systems are compared. For a dis- 
cussion of strongly intensive quantities see fl6| . 

The analysis of chemical fluctuations can be performed 
with the help of $ in two different but fully equivalent 
ways. In the first method using the identity variable, 
chemical fluctuations are treated in analogy to fluctua- 
tions of transverse momentum. In the second method <& 
is calculated from the moments of the multiplicity distri- 
butions (To| . 

We next describe the first method for the example of a 
two-component system of pions and kaons. One defines 
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the single particle variable x as x = wk, where wk is 
called the kaon identity and w K = 1 if the i— th particle 
is a kaon and w K = if the i— th particle is a pion. This 
implies unique particle identification. One then directly 
uses the definition ([3]) to evaluate $. 

Let us now discuss the most important case for which 
the measure $ of chemical fluctuations vanishes. Since 
the inclusive distribution of wk equals 



P{w K ) 



for wk = 



(N) 

W for w K = 1 



(4) 



one finds 



o 9 9 

z z = wj^ — wk 



(N k )( (N k )\ f . 



When inter-particle correlations are absent, the distri- 
bution of particle identities in events of multiplicity N 
reads 

P N (w K , O = 7VP«) F«) • • • P«) , 

(6) 

where 'P/v is an arbitrary multiplicity distribution of par- 
ticles of any type. One shows that (Z 2 ) computed with 
the event distribution © equals (Z 2 ) = (N)z 2 and con- 
sequently, $ = 0. 

In the second method the measure <£> of chemical fluc- 
tuations is obtained from the moments of the experimen- 
tally measured multiplicity distributions of kaons and pi- 
ons. As shown in Ref. [lOj . one has 



-3 (Nk)(N„) 



(AO 2 



(7) 



(A) 



((N k )-(Nk) 2 )+ { -^((N 2 )-(N^ 



(N) 3 

(A)3 



(A)s 

«JVjciV»> - (A K ){A W )) 



(8) 



which substituted in Eq. Q give the measure $. 

The formulas ([71 |SJ clearly show that $, like Vd yn , 
vanishes when the multiplicity distributions of pions and 
kaons are both poissonian and independent from each 
other. However, more generally, $ vanishes for any mul- 
tiplicity distribution provided it satisfies ([6]). The dis- 
tribution (|6|) leads to the multiplicity distribution of the 
form 



(9) 



x (Nk + NS\({Nk)\ Nk ( 1 (Nk)" N * 



(N) 



(N) 



with the moments 



V/JV, -I, =^f(A(A-l)). (ID) 



(NniN n -l)) = @£{N(N-i)) , (11) 



(AW=<^pW-l)> 



(12) 



One checks that $ and Vdyn both vanish when these mo- 
ments are substituted into Eq. (JSJ) and Eq. @, respec- 
tively. 

It appears convenient for our further considerations to 
modify <3> to the form 



* = 



(Z 2 ) 



(A) 



(13) 



which preserves the properties of $ - it vanishes in the 
absence of inter-particle correlations and it is strongly 
intensive. The measure ^ will be used to formulate the 
identity method for the study of chemical fluctuations. 
When expressed through moments of the multiplicity dis- 
tribution, it equals 



1 



(Nl){N K ) 2 + {N„) 2 (N 2 K ) 



(14) 



2{N*)(N K ){N V N K ) - (N„) 2 (N K ) - (A W )(A^) 5 



Comparing Eq. ((2} to Eq. (fl4]) one finds that ^ and fdyn 
are proportional to each other: 



*" (A)3 



(15) 



We note here that Vdyn is not intensive but it becomes 
even strongly intensive when multiplied by (A), (Nk) or 



III. EFFECT OF MISIDENTIFICATION 

As mentioned in the Introduction, complete identifi- 
cation of every particle is impossible. In this section we 
show how the incomplete particle identification influences 
the magnitudes of fluctuation measures. For this purpose 
we considered a simple model of chemical fluctuations 
where the multiplicity of pions is poissonian with a mean 
value of 100 and the number of kaons is 20% of the num- 
ber of pions (strict correlation of the numbers of kaons 
and pions) . Actually, Nk is taken as the integer number 
closest to N„/5 which is smaller or equal to N^/5. The 
fluctuation measures can be easily computed analytically 
for the model but our aim here is to simulate the effect 
of incomplete particle identification. 

There are many experimental techniques to measure 
particle mass. We discuss here the effect of misidentifi- 
cation referring to measurements of energy loss, dE/dx, 
in a detector material. This method is applied by 
the experiments NA49, NA61 and STAR. The detectors 
are equipped with Time Projection Chambers in which 
dE/dx is measured. The value of dE/dx can be used 
to identify particles because it depends both on particle 
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dE/dx [arb. units] dE/dx [arb. units] 



FIG. 1: (Color online) The distribution of dE/dx with non-overlapping peaks of pions and kaons (a), which allows unique 
particle identification, and the distribution with overlapping peaks (b), which does not allow unique identification. 



mass m and momentum p in the combination of veloc- 
ity (/3 = p/y^rri 2 + p 2 ). In the case of large separation 
of the energy loss distributions of pions and kaons, as 
schematically shown in Fig. [TJi, almost unique particle 
identification is possible. This is, however, not possible 
when the measured pion and kaon dE/dx distributions 
overlap, as illustrated in Fig. QJ). 

Performing the Monte Carlo simulations we assumed 
that the dE/dx distributions of pions and kaons are gaus- 
sians centered, respectively, at 1.4 and 1.2 in arbitrary 
units. They are normalized to the mean multiplicity of 
pions and of kaons, respectively. In order to quantify the 
bias caused by particle misidentification a simple particle 
identification scheme is used, namely, a particle is identi- 
fied as pion if dE/dx > 1.3 and as kaon if dE/dx < 1.3. 
The width a of both gaussians is chosen to be the same 
but its value is varied from to 0.08. With growing 
width of the peaks of the dE/dx distribution, the frac- 
tion of misidentified particles obviously increases. The 
results of the simulation are illustrated in Fig. [2] where 
the fluctuation measures $, ^> and i/dyn are shown as a 
function of a. As seen, the magnitudes of $, ^ and v<± yn 
decrease as the fraction of misidentified particles grows 
and the measures vanish when particle identification be- 
comes totally random. This sizable and experimentally 
unavoidable effect was the main motivation to develop 
the identity method which fully eliminates the problem. 



IV. IDENTITY METHOD 

The identity method, which is described here for a two- 
component system of pions and kaons, utilizes the mea- 
sure "J defined by Eq. (fT5|) . However, the kaon identity 



wk is not limited to either 1 or any more, but can take 
any value from the interval [0, 1]. 

In the previous section we assumed that particles 
are identified according to the energy-loss distribution. 
To make the presentation of the identity method more 
general we assume here that particle identification is 
achieved by measurement of particle mass not specify- 
ing the particular experimental technique which is used 
for this purpose. Since any measurement is of finite reso- 
lution, we deal with continuous distributions of observed 
masses of pions and kaons which are denoted as /^(m) 
and pk{iti), respectively. They are normalized as 



dmpnim) = (N n ) 



dmp K (m) = (N K ) ■ (16) 



The kaon identity is defined as: 

p K (m) 



w K (m) 



p(m) 



(17) 



where p{m) = p-nijn) + pxija), and is normalized as 

dmp(m) = (N) = + (N K ) . (18) 



If the distributions p- K (m) and p_R-(m) do not overlap, 
the particles can be uniquely identified and wk = for 
a pion and wk = 1 for a kaon. When the distributions 
p w (m) and pK(m) overlap, wk can take the value of any 
real number from [0,1]. Figure [3] illustrates the latter 
case. The mass distributions are shown in Fig. [3^, and the 
distribution of kaon identity in Fig. [3Jd. The peaks close 
to and 1 in Fig. [3h> correspond to the mass regions in 
which pions, respectively kaons, are well identified. The 
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FIG. 2: (Color online) The measures of chemical fluctuations $, \l/ (a) and Vdyn (b) as functions of the width of the energy-loss 
distribution. 



wk values around 0.5 correspond to particles for which 
the measured mass is in the transition region between the 
kaon and pion peaks (m w 280 MeV) in the distribution 
p(m) shown in Fig. [3^. 

Let us now explain how the fluctuation measure ^ 
is calculated once the mass distributions were experi- 
mentally obtained. The single particle variable entering 
Eq. (|13| is defined as in Sec. [TT) z = Wk — Wk- The 
bar denotes the inclusive average which is computed as 
follows: 



w K 



1 

JN) 
1 

(iV) 



dm p{m) WK{m) 

(Nk) 



dm pxijn) 



(N) 



Analogously one finds w K and 
quantity (Z 2 ) is obtained as 



w K 



1 Nov 



E 

i=l 



w l K - N n W K 



(19) 



The 



(20) 



where A cv is the number of events and N n is the mul- 
tiplicity of the n— th event. Substituting (Z 2 ), z 2 and 
(N) into Eq. (IT3|) . one finds the measure >F, the magni- 
tude of which, however, is biased by the effect of particle 
misidentification. Next we discuss the correction proce- 
dure. 

As shown in the Appendix, the measure \& can be ex- 
pressed through the moments of the multiplicity distri- 
butions of pions and kaons as 



* = A 



(Nk) _ 



(21) 



where 

A = 



1 

JN) 



(N 



+ (N%) - (N K ) - 



(Nk) 2 



-2(N„N K ) 



(Nk) 



and 



u K 



(N 



Kl 



dm pk (m) wk (m) 



(22) 



(23) 



In the case of complete particle identification (CI) the 
distributions p w (m) and pK^m) do not overlap and thus 
UK = 1- Then, the result for VP is: 



ci 



.4 



(Nk) 
(N) 



- 1 



(24) 



which is equivalent to the expression (1141) . 

Although particle-by-particle identification is usually 
difficult, statistical identification is reliable. In the latter 
case, we do not know whether a given particle is a kaon, 
but we know the average numbers of kaons and of pions. 
We introduce the concept of random identification which 
assumes that for every particle the probability of being a 
kaon equals (Nk)/(N). Such a situation is described by 
mass distributions of the form 



Pi(m) 





{Ni) 



m max — m n 

o 



for m < m m in , 
for m min < m < m max , 
for TO max < m , 



where i — n, K and m min and m max denote lower and 
upper limits of the measured mass range, respectively. 
With this distribution 



UK = W K 



(Nk) 
(N) 



(25) 



When u~k = (N K )/(N) is substituted into Eq. (|2H . 
"F = 0. Thus, the measure $ vanishes when particle 
identification is random. 
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FIG. 3: (Color online) The distributions of observed masses of pions p 7r (m), of kaons pK(m), and their sum p(m) (a) and the 
corresponding distribution of kaon identity wk (b) . 



Finally, we arrive at the crucial point of the consider- 
ations. It appears that the measure , F can be expressed 
as 



* = *cil 1 



YlX 
v R ) 



(26) 



where ^ci is the measure ^ for the complete identifica- 
tion, as given by Eq. (l24t . The quantities Vi and Vr are 
the values of the integral 



V— dm p(m) wk(iti) (l — w^(m)) 



(27) 



evaluated for the cases of imperfect and random identifi- 
cation, respectively. 

One proves the equality (|2"?J)) by observing that 



Vi = {N K ){i-m) 



Vr = 



(Nk)(N v ) 
(N) 



(28) 



Substituting Eqs. (JHHS]) into the equality (J2H]), one ob- 
tains the formula (|21l) . Actually, the relation (f2l)|) was 
first discovered by performing various numerical simula- 
tions and only then it was proven analytically. Equa- 
tion (|26[) allows one to experimentally obtain <J/ci from 
fy. Thus, the effect of misidentification is fully corrected 
by the factor (1 - Vj/Vr) 2 which measures the quality 
of the applied procedure of particle identification. The 
factor is independent of the correlations under study and 
can be determined from experimental data. 



procedure to obtain the measure Vfci of fluctuations of 
kaons (or any other selected particle type) with respect to 
all particles (represented by the sum of kaons and pions 
in the previous sections). It is important to note that if 
misidentification occurs between more than two particle 
types, the identity method does not allow to study rela- 
tive fluctuations of two of them, e.g. of kaons and pions 
in the presence of protons. 

We come back to the specific, but typical, example 
considered in Sec. Ill of particle identification via mea- 
surements of particle energy loss in the detector material. 
The energy loss dE/dx is denoted by X. The energy 
loss of particles of a given type depends on the parti- 
cle mass and momentum (via the velocity) and detector 
characteristics. The distribution of X is usually a multi- 
dimensional function, which can be determined experi- 
mentally by averaging over particles from many interac- 
tions. The energy loss distribution, which for a given 
particle momentum is typically fitted by a sum of four 
gaussians corresponding to electrons, pions, kaons and 
(anti-)protons, allows one to determine the average mul- 
tiplicities of kaons and of all particles. Having obtained 
this information, one should proceed as follows. 

(i) Extract the energy- loss distribution of kaons pk (X) 
from the inclusive distribution p(X). The distribu- 
tions should be normalized as 



dX p(X) = (N) , 



dX PK {X) = {N K ) 



EXPERIMENTAL PROCEDURE 



(ii) Determine the kaon identity 



The application of the identity method to experimen- 
tal data is not difficult. We present here a step-by-step 
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for every registered particle. 

(iii) Compute the fluctuation measure "J from the defi- 
nition (|13[) . This is the raw value of which is not 
corrected yet for the effect of misidentification. 

(iv) Knowing the mean multiplicities, compute the 
quantity 



Vr = 



(N K )((N) - (N K )) 
(N) 



(v) Using all particles calculate the integral 



Vi 



dXp{X)w K {X)(l-w K {X)) 



(vi) Determine the fluctuation measure ^cij which is 
free of the effect of misidentification, as 



(1 - ViIVrY 



The experimental data accumulated by the NA49 Col- 
laboration are currently under analysis using the identity 
method proposed in this paper (l7j . 
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Appendix 



We express here the measure ^> through the moments of multiplicity distributions of pions and kaons. For this 
purpose one writes (Z 2 ) as 



OO OO p p p 



(A.l) 



N„=0 N K =0 



x / dmfPjf(mf) / dmfP K {m«)--- I dm* K P K {m% K ) 



.K 



.K 



x (w K (mJ) + w K {ml) + ...+ w K (mJfJ + w K (m^) + w K {m%) + ... + w K (m§ K ) - (N* + N K )w K ^j , 

where Vn^Nk IS the multiplicity distribution of pions and kaons; P^im) = p n (m) / '{N n ) and Pr-(to) = pk(iti)/ (Nk) 
are the mass distributions of pions and kaons, respectively. 
Equation (|A.ip gives 



(Z 2 ) = (N n )ul + (N k )u% + (N 2 }wk 2 + (N„{N„-1))T^ 2 + (N k {N k -1))uk 2 
+ 2(N^N K )u^uW - 2(NN- K )w K ~u^ ~ 2{NN k )wFuk , 



(A.2) 



where 



(Nk) 



dm pk(iti) w K (m) , 



with n = 1,2. Since 



(NJ (N K )— — 



Equation (|A.2[) provides 



(AO 5 



(Z 2 ) = (N)w 2 K + [(N 2 ) + (N w (N w -l))-^-2(NN«) 



(N) 



w K 



(A.3) 



(A.4) 



(A.5) 
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- 2 



(n A n, - u> W££> _ w> w + {NNk) _ {NM sm 



w K u K 



(NAN, - 1)>^2- + (Nk(N k - 1)) - 2(N,N K ) {Nk) 



U K 



Keeping in mind that z 2 = w k — wk 2 and u>k — (Nk)/(N), one finds after somewhat lengthy calculations the formula 
(f2"Tj) with A given by Eq. 
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